^I^tty. Docket No. CLG 99-002 

APPARATUS AND METHODS FOR MEASUREMENT OF DENSITY OF MATERIALS 
USING A NEUTRON SOURCE AND TWO SPECTROMETERS 

BACKGROUND OF THE INVENTION 

FIELD OF THE INVENTION 

This invention is directed toward measurement of density of material, and more 
particularly directed toward a system for measuring bulk density of material penetrated by a 
borehole, wherein the system comprises a source of neutron radiation and preferably two radiation 
detection spectrometers. The system can alternately be embodied to measure other material 
properties and to measure density of materials not penetrated by a borehole. 

BACKGROUND OF THE ART 

Systems utilizing a source of radiation and a radiation detector have been used in the prior 
art for many years to measure density of material. One class of prior art density measuring 
systems is commonly referred to as "transmission" systems. A source of nuclear radiation is 
positioned on one side of material whose density is to be measured, and a detector which responds 
to the radiation is positioned on the opposite side. After appropriate system calibration, the 
intensity of measured radiation can be related to the bulk density of material intervening between 
the source and the detector. A second class of prior art density measuring systems is commonly 
referred to as "back scatter" systems. Both a source of nuclear radiation and a detector, which - 
responds to the radiation, are positioned on a common side of material whose density is to be 
measured. Radiation impinges upon and interacts with the material, and a portion of the 
impinging radiation is scattered by the material and back into the detector. After appropriate 
system calibration, the intensity of detected scattered radiation can be related to the bulk density of 
the material. 

Backscatter type systems have been used for decades to measure density of material, such 
as earth formation, penetrated by a borehole. The measuring instrument or "tool" typically 
comprises a source of radiation and at least one radiation detector which is axially aligned with the 
source and typically mounted within a pressure tight container. 
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Systems which employ the backscatter configuration with a source of gamma radiation and 
one or more gamma ray detectors are commonly referred to as "gamma-gamma" systems. Sources 
of gamma radiation are typically isotopic such as cesium- 137 ('^^Cs), which emits gamma 
radiation with energy of 0.66 million electron volts (MeV) with a half life of 30.17 years. 
Alternately, cobalt-60 (^^Co) is used as a source of 1 . 1 1 and 1 .33 MeV gamma radiation with a 
half life of 5.27 years. The one or more gamma ray detectors can comprise ionization type 
detectors, or alternately scintillation type detectors if greater detector efficiency and delineation of 
the energy of measured scattered gamma radiation is desired. 

The basic operational principles of prior art gamma-gamma type backscatter density 
measurement systems are summarized in the following paragraph. For purposes of discussion, it 
will be assumed that the system is embodied to measure the bulk density of material penetrated by 
a borehole, which is commonly referred to as a density logging system. It should be understood, 
however, that other backscatter density sensitive systems are known in the prior art. These 
systems include tools which use other types of radiation sources such as neutron sources, and 
other types of radiation detectors such as detectors which respond to neutron radiation or a 
combination of gamma radiation and neutron radiation. 

A backscatter gamma-gamma density logging tool is conveyed along a well borehole 
penetrating typically earth formation. Gamma radiation fi-om the source impinges upon material 
surrounding the borehole. This gamma radiation collides with electrons within the earth formation 
material and loses energy by means of several types of reaction. The most pertinent reaction in 
density measurement is the Compton scatter reaction. After undergoing typically multiple 
Compton scatters, a portion of the emitted gamma radiation is scattered back into the tool and 
detected by the gamma radiation detector. The number of Compton scatter collisions is a fimction 
of the electron density of the scattering material. Stated another way, the tool responds to electron 
density of the scattering earth formation material. Bulk density rather than electron density is 
usually the parameter of interest. Bulk density and electron density are related as 

(1) Pe = Pb(2(IZi)/MW) 




where 
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Pe = the electron density index; 
Pb = the bulk density; 

(SZi) = the sum of atomic numbers Zi of each element i in a molecule of the material; and 
MW = the molecular weight of the molecule of the material. 

For most materials within earth formations, the term (2 (SZi) / MW) is approximately equal to 
one. Therefore, electron density index pe to which the tool responds can be related to bulk density 
Pb, which is typically the parameter of interest, through the relationship 

(2) Pb = Ape + B 

where A and B are measured tool calibration constants. Equation (2) is a relation that accounts for 
the near linear (and small) change in average Z/A that occurs as material water fraction changes 
with material porosity, and hence changes with bulk density. 

The radial sensitivity of the density measuring system is affected by several factors such as 
the energy of gamma radiation emitted by the source, the axial spacing between the source and 
one or more gamma ray detectors, and properties of the borehole and the formation. Formation in 
the immediate vicinity of the borehole is usually perturbed by the drilling process, and more 
specifically by drilling fluid "invades" the formation in the near borehole region. Furthermore, 
particulates from the drilling fluid tend to buildup on the borehole wall. This buildup is 
commonly referred to as "mudcake". Mudcake, invaded formation and other factors perturbing 
the near borehole region can adversely affect a formation bulk density measurement. It is of prime 
importance to maximize the radial depth of investigation of the tool in order to minimize the 
adverse effects of near borehole conditions. Generally speaking, an increase in axial spacing 
between the source and the one or more detectors increases radial depth of investigation. 
Increasing source to detector spacing, however, requires an increase in source intensity in order to 
maintain acceptable statistical precision of the measurement. Prior art systems also use multiple 
axial spaced detectors, and combine the responses of the detectors to "cancel" effects of the near 
borehole region. This method is marginally successfiil since nuclear systems are inherently 
shallow depth of investigation. Depth of investigation can be increased significantly by increasing 
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the energy of the gamma-ray source. This permits deeper radial transport of gamma radiation into 
the formation. Unfortunately, there are no isotopic sources emitting gamma radiation above 1.33 
MeV which have a half-life sufficiently long for typical commercial use and which are reasonably 
inexpensive to produce. Accelerator sources have been used in the prior art to generate gamma 
radiation of energy greater than 10 MeV. These sources are, however, physically large, costly to 
fabricate, costly to maintain, and often not suited for harsh environments such as a well borehole. 

SUMMARY OF THE INVENTION 

This invention is directed toward a system for measuring density and other properties of 
material penetrated by a borehole. Alternately, the system can be embodied for material analysis 
in a variety of non-borehole envirormients. Configuration of the system is based upon the 
backscatter concept discussed in the previous section of this disclosure. 

The sensor instrument or "tool" comprises preferably an axially spaced source of radiation 
and preferably two axially spaced radiation detectors, which discriminate energy of radiation 
impinging upon the detectors. 

The source is preferably a neutron source which emits, or induces within material being 
measured, gamma radiation with energy greater than energy obtainable with isotopic gamma ray 
sources. Further, the neutron-induced gamma source is dispersed within the material being 
measured, providing for a larger investigation depth than obtainable with prior art systems 
wherein a gamma ray source is located within a tool and not dispersed within material being 
measured. The dispersed source of gamma rays results fi-om the reaction of transported neutrons 
with the materials surrounding the tool. The neutron source is preferably isotopic, although other 
types of radiation sources such as a neutron generator can be used. Conceptually, a source of 
gamma radiation can be used but an accelerator type gamma ray source would be required to 
obtain the desired high energy radiation. Furthermore, such a source has the disadvantage of 
being a point source outside of the material being measured rather than a dispersed gamma source 
within the material being measured. 

A first suitable isotopic neutron source is americium 241 -beryllium (AmBe), which 
produces alpha, gamma and neutron radiation. More specifically, the AmBe source produces 
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gamma radiation at 4.43 MeV from the decay of carbon- 12 ('^C), which is produced by alpha 
radiation from americium interacting with beryllium. Alternately, plutonium-beryllium (PuBe) 
can be used with similar results. Details of the alpha-beryllium reaction will be discussed in a 
subsequent section of this disclosure. The 4.43 MeV carbon gamma radiation is much more 
energetic that previously discussed gamma radiation from '^^Cs (0.66 MeV) and ^°Co (1.11 and 
1.33 MeV). A gamma-gamma back scatter density measuring system using 4.43 MeV radiation 
will obtain much greater penetration into material being measured that that obtained with the prior 
art '^'Cs or ^Co sources. More importantly, however, AmBe induces a source of gamma 
radiation within the material being measured resulting in a deeper depth of investigation. 
Neutrons from the AmBe source enters the material being measured. The neufron radiation 
produces gamma radiation primarily by inelastic scatter and thermal capture reactions with nuclei 
within the material. This neufron induced gamma radiation, and to a lesser proportion the back 
scattered 4.43 MeV gamma radiation, is detected by the preferably two axially spaced scintillation 
type gamma radiation specfrometers. The specfrometer positioned closest to the source is 
hereafter referred to as the "short spaced" detector, and the specfrometer positioned farthest from 
the source is hereafter referred to as the "long spaced" detector. Other types of detectors 
exhibiting spectral gamma ray response can be used. Alternate detectors include solid state 
detectors and gas filled detectors fabricated to be energy dependent. 

A second suitable isotopic neufron source is Califomium-252 (^^^Cf). ^^^Cf does not 
produce 4.43 MeV gamma radiation internally as does AmBe (and PuBe). Neufrons are, however, 
produced, enter the material being measured, and produce an (n-y) radiation with energy specfral 
characteristics similar to the (n-y) component generated by AmBe. Energies of a significant 
portion of the induced (n-y) radiation are greater than energies obtainable with isotopic gamma ray 
sources. The induced gamma radiation is detected by the short spaced and the long spaced 
detectors. Fission gamma radiation is also produced and incorporated into the measurement. 

The system embodied as a borehole device is conveyed along a borehole by means of a 
wireline or a drill string. Gamma ray specfra are measured in both long spaced and short spaced 
detectors. Spectra from each detector are split and integrated over preferably two energy regions 
or energy "windows". The energy windows are preferably contiguous, with a low energy window 
extending from a few hundred keV to about 3.0 MeV, and the high energy window extending 
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from about 3.0 MeV to about 7 to 10 MeV. Energy windows for the short spaced and long spaced 
detectors are preferably identical, although different window limits can be employed using 
appropriate normalization methods. It should be understood that the energy limits of both the low 
window and the high window can be varied as long as the low energy window is responsive to 
down scatter gamma radiation and the high energy window is relatively insensitive to down 
scattered gamma radiation. As an example, a low energy window extending from several hundred 
keV to about 2 MeV would encompass a large percentage of down-scattered gamma radiation and 
backscatter radiation. 

A normalization factor is next computed, wherein counts recorded in the high energy 
window of the long spaced detector multiplied by the normalization factor equals the counts 
recorded in the high energy window of the short spaced detector. Counts recorded in the low 
energy window of the long spaced detector are multiplied by the normalization factor thereby 
yielding a normalized low energy counts for the long spaced detector. The count measured in the 
low energy window of the short spaced detector is then subtracted from the normalized low energy 
count from the long spaced detector thereby yielding a low window count difference. For a given 
bulk density of material being measured, the percent of down scatter radiation increases with 
detector spacing. Furthermore, the percent of down scatter radiation is proportional the bulk 
density of the material being measured. These factors combine so that the low window count 
difference varies as a function of material bulk density, with low energy count difference 
increasing with increasing material bulk density. This process is recited mathematically in 
equations (1) and (2) discussed previously. 

The fiinctional relationship between low window count difference and material bulk 
density is determined by (1) calibrating by measuring tool response in materials of known density, 
or (2) by suitable tool response calculations such as Monte Carlo calculations, or (3) by a 
combination of measured and calculated tool responses. 

BRIEF DESCRIPTION OF THE DRAWINGS 



So that the manner in which the above recited features, advantages and objects of the 
present invention are obtained and can be understood in detail, more particular description of the 
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invention, briefly summarized above, may be had by reference to the embodiments thereof which 
are illustrated in the appended drawings. 

Fig. I illustrates the invention embodied as a wireline logging system; 
Fig. 2 illustrates the invention embodied as a logging-while-drilling system; 

Fig. 3 a is a gamma ray spectrum measured with the short spaced detector using Cf as a 
source of radiation; 

Fig 3b is a gamma ray spectrum measured with the long spaced detector using ^^^Cf as a 
source of radiation; 

Fig. 3c is a gamma ray spectrum measured with the short spaced detector using AmBe as a 
source of radiation; 

Fig. 3d is a gamma ray spectrum measured with the long spaced detector using AmBe as a 
source of radiation; 

Fig. 4a is a plot of low window count difference as a function of gamma ray energy for 
differing density material and using ^^^Cf as a source of radiation; 

Fig. 4b is a plot of low window count difference as a function of gamma ray energy for 
differing density material and using AmBe as a source of radiation; 

Fig. 5 is a conceptual representation of the functional relationship between low window 
count difference and bulk density for ^^^Cf and AmBe sources; 

Fig. 6 is a flow chart of the data processing method; and 

Fig. 7 depicts the invention embodied to measure material bulk density and other material 
properties in a non-borehole geometry. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

This disclosure sets forth a system comprising apparatus and methods for measuring 
density and other properties of material. Other properties include material porosity and lithology. 

The system can be embodied to measure properties of earth formation penetrated by a 
borehole. A measuring instrument or "tool" is conveyed within the borehole thereby measuring 
parameters of formation penetrated by the borehole. Properties are typically measured as a 
function of tool depth within the borehole thereby generating a "log" of the measured parameter. 
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The tool can be conveyed by means of a wireline or a drill string. Alternately, the system can be 
embodied for material analysis in a variety of non-borehole environments. Configuration of the 
system is based upon the measure of back-scattered nuclear radiation as discussed conceptually in 
a previous section of this disclosure. Limitations of prior art backscatter systems are minimized 
using the present invention, as will become apparent in subsequent sections of this disclosure. 
Depth of investigation of the measurement is improved substantially by the dispersed nature of the 
gamma source as produced by neutron reactions in the surrounding materials. Investigation depth 
is also increased by the use of gamma radiation of energy higher than that available fi-om common 
isotopic gamma-ray sources. 

Fig. 1 illustrates the invention embodied as a wireline well logging system. The system 
comprises a logging tool 10 connected to a wireline logging cable 36 by means of a connector 34. 
The tool 10 is shown suspended within a borehole 13 defined by a borehole wall 12 and 
penetrating an earth formation 14. The logging cable 36 passes over a sheave wheel 40 at the 
surface of the earth 43. A depth measuring means 42 cooperates with the sheave wheel 40 to 
generate a signal indicative of the depth of the tool 10 within the well borehole 13. Response 
signals from the tool 10 are processed by a first spectrometer processor 46 and a second 
spectrometer processor 44. Outputs from the first and second spectrometer processor 46 and 44 
are input along with a signal fi-om the depth measuring means 42 into a central processing unit 
(CPU) 48. Formation bulk density and other formation and borehole parameters of interest are 
computed within the CPU 48 and output to a recording means 50 which can comprise a digital 
memory, a digital recording means, or hard copy record as identified by the numeral 52. 

Attention is now directed to the various elements of the tool 10 as shown in Fig. 1. The 
tool 10 comprises a source of radiation 20. The source 20 is preferably an isotopic neutron source 
such as Califomium-252 (^^^Cf), or alternately americium-beryllium (AmBe) which produces 
alpha, gamma and neutron radiation. The source 20 is surrounded by a fast neutron moderating 
material 18 such as high density polyethylene (HDPE) to maximize the thermal neutron flux 
entering the formation 14. 

12 

Americium emits alpha particles which react with Be to form C through the reaction 
(3) ^Be(a,n)'^C* . 
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thereby producing neutrons over a wide energy range. Excited C* decays to the ground state by 
the emission of 4.43 MeV gamma radiation. 

^^^Cf produces fast neutrons by spontaneous fission. Gamma rays are also produced by the 
spontaneous fission of ^^^Cf The average neutron energy fi-om ^^^Cf is about 2 MeV. The 
average gamma energy is about 0.7 MeV. The average energy of neutrons produced by an AmBe 
source from the Am (a,n)Be reaction is about 4.5 MeV. In addition to differences in average 
energy, the neutron energy distribution of the two sources is significantly different. ^^^Cf neutrons 
are broadly peaked at about 1 MeV. Above about 2 MeV, neutron intensity decreases 
exponentially with increasing energy by about a factor of 10 every 3 MeV. On the other hand, 
AmBe neutrons exhibit several energy peaks in the range from a few hundred kilo elecfron volts 
(keV) to about 10 MeV. This energy distribution difference for the Cf and AmBe neutron 
sources results in differing formation neutron distributions for a given formation and borehole 
condition. These differing neutron distributions result in differing neutron induced gamma ray 
spectra. The fission gamma radiation is generally low in average energy and decreases in intensity 
with increasing energy. This fission gamma radiation can be effectively biased out by the setting 
of the lower limits of the low energy windows of the detectors. 

A short spaced detector comprising a scintillation crystal 24 and a photomultiplier tube 26 
is positioned a distance 30 from the source 20, and shielded from the source with a gamma ray and 
neufron shielding material 16 such as tungsten gammaloy. The shielding 16 minimizes neufron 
and gamma radiation from "streaming" directly from the source 20 into the scintillation crystal 24. 
Fast neutron moderating material 1 8 surrounds the source 20 to increase the flux of thermalized 
neufrons entering the formation 14. The scintillation crystal 24 is preferably sodium iodide (Nal), 
but other scintillation crystal materials can be used. The scintillation crystal 24 is radially shielded 
with a material 22 which readily absorbs thermal neutrons impinging thereon thereby preventing 
thermal neutron activation and capture reactions from occurring within the scintillation crystal. A 
suitable material for the shielding material 22 is HDPE, which moderates fast neufrons, and 
including additives such as Cd, Gd and B which are all prolific thermal neufron absorbing 
elements. 

Still referring to Fig. 1, a long spaced detector comprising a scintillation crystal 25 and 
photomultiplier tube 27 and which is axially spaced from the source 20 by a distance 3 1 . Neutron 
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and gamma ray shielding 16 are disposed between short spaced and the long spaced detectors. 
The spacings 30 and 31 are preferably about 7 inches and 14 inches, respectively. These spacings 
are based on laboratory measurements. The long spacing is preferably greater than 14 inches for 
better density sensitivity, but will be limited by source strength and corresponding count rate 
statistics considerations. The shielding-spacing arrangement minimizes response of the short and 
long spaced detectors to direct flux from the source 20, and also maximizes the desired signal 
related to formation bulk density as will be discussed in a subsequent section of this disclosure. It 
is preferred that the long and short spaced detectors be identical in response characteristics to 
facilitate normalization, which will be discussed in a subsequent section of this disclosure. 
Alternately, different long and short spaced detector responses can be incorporated using 
appropriate response normalization methods. All elements of the tool 10 are mounted within a 
pressure housing 1 7, and powered and controlled by an electronics package 32. 

Fig. 2 illustrates conceptually the tool 10 embodied as a logging-while-drilling (LWD) 
system. The tool 10 is conveyed by a drill string 36' along the borehole 13 penetrating the 
formation 14. A drill bit 11 terminates the drill string 36'. Surface drilling and recording 
apparatus are well known in the LWD art, and are therefore illustrated conceptually at 41 for 
brevity. 

Embodied either as a wireline or LWD system, neutrons and gamma radiation from the 
source 20 enter the formation 14. The neufron radiation produces gamma radiation by inelastic 
scatter and thermal capture reactions with nuclei within the formation 14. A portion of the neufron 
induced radiation is down scattered as it is fransported through the formation 14 and impinges 
upon the long and short spaced detectors. 

Fig. 3a is a spectrum 62 of gamma ray intensity plotted as a fiinction of gamma ray energy 
as measured by the short spaced detector using Cf as a source 20. The peak 60 at 2.23 MeV is 
produced by thermal capture of neufrons by hydrogen. The peak at 63 is produced by inelastic 
scatter of fast neutrons by the element silicon (Si). Measures of counts in these peaks can be 
related to other parameters of the material being measured as will be discussed in subsequent 
sections of this disclosure. These and other peaks are superimposed on a substantial and complex 
distribution of inelastic scatter and capture gamma radiation induced by neutrons from the source 
20 interacting with nuclei within the formation 14. The spectrum is partitioned into two 
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continuous low and high energy windows. A low energy window extends from gamma ray energy 
of a few hundred keV to an energy of about 3.0 MeV, and is designated by the numeral 68. An 
altemate upper energy of about 2 MeV (not shown) can be used to exclude any unscattered 2.223 
MeV hydrogen thermal capture gamma radiation. A high energy window extends from about 3.0 
MeV to an energy of about 7.0 MeV or greater. Counts recorded within the low energy window 
are identified by the shaded area 71 and are hereafter designated as Cl,ss, where the "L" identifies 
the low energy window and "SS" identifies the short spaced detector. Counts recorded within the 
high energy window are identified by the shaded area 70 and hereafter designated as Ch,ss, where 
the "H" identifies the high energy window and "SS" again identifies the short spaced detector. 
Counts in the peaks 60 and 63 can be obtained by measuring counts in the energy windows 60' and 
63', respectively, or by other methods such as spectrum stripping. 

Fig. 3b is a spectrum 62' of gamma ray intensity plotted as a fianction of gamma ray energy 
as measured by the long spaced detector using ^^^Cf as a source 20. The peak 60 again identifies 
the 2.223 MeV gamma ray produced by thermal capture of neutrons by hydrogen, which again is 
superimposed on a substantial and complex distribution of inelastic scatter and capture gamma 
radiation induced by neutrons from the source 20 interacting with nuclei within the formation 14. 
The spectrum 62' is lower in intensity than the short spaced detector spectrum 62 shown in Fig. 3 a 
due to a greater source-detector spacing 3 1 . The long spaced detector spectrum 62' is partitioned 
once again into two continuous low and high energy windows extending from near a few hundred 
keV to 3.0 MeV and 3.0 MeV to about 7.0 MeV or larger, respectively. Counts recorded within 
the low energy window are identified by the shaded area 73 and are hereafter designated as Cl,ls5 
where the "L" identifies the low energy window and "LS" identifies the long spaced detector. 
Counts recorded within the high energy window are identified by the shaded area 72 and are 
hereafter designated as Ch,ls, where the "H" identifies the high energy window and "LS" again 
identifies the long spaced detector. 

A normalization factor is next computed from the relationship 

(4) N = Ch,ss/Ch,ls 

where Ch,ss and Ch,ls have been defined previously. Physically, equation (4) normaUzes the long 
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spaced and short spaced detectors to counts recorded in their high energy windows. Counts 
recorded in the low energy window of the long spaced detector, Cl,ls, are multiplied by the 
normalization factor N thereby yielding a normalized low energy count for the long spaced 
detector which is 



The count measured in the low energy window of the short spaced detector, Cl,ss, is then 
subtracted from the normalized low energy count for the long spaced detector, C'l,lsj thereby 
yielding a low window count difference 



For a given bulk density of material being measured, the percent of down scatter radiation 
increases with detector spacing. Stated another way, percent down scatter is greater in the long 
spaced detector. Furthermore, the percent of down scatter radiation is proportional the bulk 
density of the material being measured. Both the proportional relationship between Compton 
scattering and electron density, and the near proportionality of electron density and bulk density, 
have been discussed in detail in previous sections of this disclosure. These factors combine so that 
the low window count difference varies as a function of material bulk density, with low energy 
count difference increasing with increasing material bulk density. This effect is illustrated 
graphically in Fig. 4a, which is a plot of low window count difference, AC, as a function of 
gamma ray energy increments or energy "channels", with ^^^Cf as a source 20. Curve 76 is a 
spectrum measured in a silicate formation (gabbro) with bulk density pb = 3.01 grams per cubic 
centimeter (gm/cm''). Curve 78 is a spectrum measured in a silicate formation (sandstone) with 
bulk density pb = 2.01 gm/cm^. Because of the normalization procedure discussed above, curves 
76 and 78 above the energy 66 (high energy window) essentially overlay at a value of AC = 0 as 
identified at 80. Below the energy 66 (low energy window), the curves deviate with curve 76 
representing the more dense material being larger. 



(5) 



C'l,ls - N Cl,ls 



(6) 



AC - Cl,ls - Cl,ss 
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The functional relationship f(Ki, AC) between low window count difference AC and 
material bulk density pb is defined as 

(7) pb = f(Ki,AC) 

where AC is the sum of values AC over a given energy window and where Kj are tool calibration 
constants determined by (1) measuring tool response in materials of known density, or (2) by 
suitable tool response calculations such as Monte Carlo calculations, or (3) by a combination of 
measured and calculated tool responses. The relationship between the quantity of interest pb and 
the quantity AC computed from measured quantities is shown conceptually in Fig 5, with f(Ki, AC) 
being represented by a curve 82 when ^^^Cf is used as a neutron source 20. 

As mentioned previously, AmBe is an alternative source 20. Fig. 3c is a spectrum 61 of 
gamma ray intensity plotted as a function of gamma ray energy for gamma radiation measured by 
the short spaced detector. AmBe is used as a source 20. The peak 60 at 2.23 MeV is again 
produced by thermal capture of neutrons by hydrogen, the peak 63 at 1 .78 MeV is produced by 
inelastic scatter of neutrons from Si, and the peak 64 at 0.85 MeV is produced by the inelastic 
scatter of neutrons fi^om iron (Fe). All peaks are superimposed on a substantial and complex 
distribution of inelastic scatter and capture gamma radiation induced by neutrons fi-om the AmBe 
source 20 interacting with nuclei within the formation 14. As discussed in the section using ^^^Cf 
as a source, the spectrum is partitioned into two contiguous low and high energy windows with 
energy ranges defined previously. Counts recorded within the low energy window are identified 
by the shaded area 71' and again designated as Cl,ss- Counts recorded within the high energy 
window are identified by the shaded area 70' and again designated as Ch.ss- 

Figs. 3d is a spectrum 61' of gamma ray intensity plotted as a function of gamma ray 
energy for gamma radiation measured by the long spaced detector using AmBe as a source 20. 
The peak structure is the same as that in the short spaced detector spectrum 61 shown in Fig. 3c, 
but the long spaced detector spectrum 61' is again lower in intensity than the short spaced detector 
spectrum due to a greater source-detector spacing 3 1 . Previously discussed peaks 60 and 64 fi-om 
Fe and Si are somewhat less pronounced. The long spaced detector spectrum 61' is partitioned 
once again into two contiguous low and high energy windows. Counts recorded within the low 
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energy window are identified by the shaded area 73' and are again designated as Cl,ls- Counts 
recorded within the high energy window are identified by the shaded area 72' and again designated 
as Ch,ls- 

The quantities N, C'l,ls and AC using AmBe as a source 20 are again computed using 
equations (4), (5) and (6), respectively. Fig. 4b is a plot of low window count difference, AC, as a 
function of gamma ray energy "channels", with AmBe as a source 20. Curve 76' is a spectrum 
measured in the silicate formation with bulk density pb = 3.01 gm/cm^. Curve 78' is a spectrum 
measured in the silicate formation with bulk density pb = 2.01 gm/cm"'. As in Fig. 4a, the 
normalization procedure is such that curves 76' and 78' above the energy 66 (high energy window) 
essentially overlay at a value AC - 0 as identified at 80. Below the energy 66 (low energy 
window), the curves deviate with curve 76' representing the more dense material being larger. It 
should be noted, however, that deviation or "spread" between curves 76' and 78' is not as great as 
deviation between curves 76 and 78 shown in Fig. 4a. The fiinctional relationship f(Ki, AC) 
between pb and AC, for AmBe as a source, is shown conceptually as curve 82' in Fig. 5. As 
previously defined, AC is the sum of values AC over a specific energy window. Physically, the 
values of Ki are such that f(Ki, AC) is not as sensitive to AC as a function of pb when AmBe is 
used as a source. This indicates that, for a given tool geometry, ^^^Cf is the preferred source 20 for 
the tool 10. 

The preference of a ^^^Cf source is further illustrated in the following tables. Tabulated 
counts are measured in known test formations and are tabulated in count rates in counts per 
second. Table 1 is a tabulation of all measured and computed parameters discussed above in 
formations of 2.01 gm/cm^ and 3.01 gm/cm^ using ^^^Cf as a source. Note that the difference or 
"spread" in AC between the two formations is AC3.01 - AC2.01 = 1397 counts per second with an 
average deviation of 281 counts per second. 
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TABLE 1 

Source: ^"Cf 

N_ fib Cl.ss C l.ls C h.ss C h.ls Q.lxs 
4.999 2.01 6588 1447 614 123 7238 650 
6.980 3.01 7323 1342 698 100 698 2046 
AC3.oi-AC2.oi= 1397 
Average deviation = 281 
Percent deviation = 20% 



10 D Table 2 is a tabulation of all measured and computed parameters discussed above in formations of 
■g 2.01 gm/cm^ and 3.01 gm/cm^ using AmBe as a source. Note that the difference or "spread" in 
Th AC between the two formations is AC3.01 - AC2.01 = 735 counts per second with an average 
deviation of 362 counts per second. 

15 L TABLE 2 

^ Source: AmBe 

iO. I^L fib Cl.ss C l,ls C h.ss C h.ls C^lxs AC 

y 5.652 2.01 12766 2417 1254 222 13659 894 

7.184 3.01 14859 2295 1454 202 16487 1629 
20 AC3.01 - AC2.01 = 735 

Average deviation = 362 
Percent deviation = 49% 



In comparing data listed in Tables 1 and 2, it is apparent that the sensitivity of the bulk density 
25 measurement, the average deviation, and a corresponding percent average deviation, are all 

superior using Cf as a source. Stated another way, the ^"Cf source yields a more sensitive 

measurement with less statistical uncertainty. As discussed previously, the superiority of the Cf 

source can also be seen graphically in comparing Figs. 4a and 4b. 

Calibration curves 82 and 82' depicted in Fig. 5 are typically for a specific lithology such 
30 as silicates or calcites. In order to maximize the accuracy of the formation bulk density 
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measurement, some knowledge of formation lithology is required. This information can be 
obtained from additional types of measurements using, as examples, acoustic or electromagnetic 
systems. Lithology information can also be extracted from the response of the present invention. 
Thermal capture or inelastic scatter gamma radiation from identifying elements of formation 
lithology can usually be identified and measured in the spectra recorded by the short spaced and 
long spaced detectors. As an example, characteristic radiation from elemental silicon as shown in 
the peak 63 of Figs. 3a and 3c can be used to identify the material being measured as containing 
silicates, and to quantify the amount of silicate from counts recorded in the energy window 63'. A 
measure of counts in one or more peaks from calcium (not shown) can be used to identify the 
lithology of the material being measured as a calcite. Additional energy windows, or other 
techniques such as spectrum fitting, can be used to quantify this characteristic radiation which is 
subsequently used to determine lithology. Corrections for formation lithology can be made 
directly in the data processing methodology. 

It is noted that the amount of hydrogen contained within the material being measured can 
be obtained from the magnitude of counts in the peak 60 using the energy window 60'. If other 
material parameters are known, such as the type of liquid saturating pores pace of the material, the 
porosity of the material can be obtained from a measure of hydrogen content. 

The data processing methodology is summarized in the flow chart of Fig 6. The process 
works equally well for data measured using AmBe or ^^^Cf as a neutron source. Parameters Cl,ss, 
Cl.ls, Ch,ss, and Ch,ls are measured at steps 84, 88, 86 and 90, respectively. The normalization 
factor N is computed at step 92. The normalized low energy window count for the long spaced 
detector, C'l,ls, is computed at step 94. The parameter AC is computed at step 96. An optional 
lithology correction, using elemental concentration measurements of the material, is made at step 
97, and AC is converted to bulk density pb at step 98. 

The CPU 48 illustrated in Fig. 1 is programmed to perform computations to convert 
measured quantities Cl,ss, Cl,ls, Ch,ss, and Ch.ls, and optionally lithology measurements and 
corrections, into one or more parameters of interest such as bulk density pb. Computations are 
summarized in the flow chart of Fig. 6. It should be understood by those of ordinary skill in the 
art that the mathematical algorithms illustrated and discussed above represent one method for the 
conversion of measured or computed parameters into one or more parameters of interest. Other 
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algorithms and methodology can be developed to obtain similar results. As an example, other 
normalization techniques can be used to normalize long and short spaced detector responses if 
high and low energy windows differ in energy range. Furthermore, the high and low windows can 
be separated or even overlapped. Differing long and short spaced detector sensitivities will also 
require modifications in the normalization process. The response of the system to changes in 
density can be reflected in other combinations of data, such as the ratio C'l,ls / Cl,ss instead of the 
difference C'l,ls - Cl,ss- The functional relationship between measured parameters, computed 
parameters, and formation parameters of interest can, however, be effectively determined using 
algorithms specifically set forth in this disclosure. 

The invention can be embodied with a non-borehole tool 10' as shown in Fig. 7. The tool 
10' comprises a source of radiation 120, which is displaced radially fi-om the major axis of the tool. 
The source 120 is again preferably ^^^Cf and is surrounded by a fast neutron moderating material 
118. A short spaced detector comprising a scintillation crystal 124 and a photomultiplier tube 126 
is axially spaced at a distance 130 fi^om the source 120, and is shielded with a gamma ray and 
neutron shielding material 116 such as tungsten gammaloy. As in the borehole embodiment, the 
shielding 116 minimizes neutron and gamma radiation fi-om "streaming" directly fi-om the source 
120 into the scintillation crystal 124. The fast neutron moderating material 118 surrounds the 
source 120 to increase the fiux of thermalized neutrons entering the material 1 14 to be measured. 
The scintillation crystal is preferably sodium iodide (Nal), but other scintillation crystal materials 
can be used. The scintillation crystal 124 is radially shielded with a material 122 which readily 
absorbs thermal neutrons impinging thereon thereby preventing thermal neutron activation and 
capture reactions from occurring within the scintillation crystal 124. A suitable material for the 
shielding material 122 is again HDPE, which moderates fast neutrons plus additives such as Cd, 
Gd and B which are all prolific thermal neutron absorbing elements. 

Still referring to Fig. 7, a long spaced detector comprising a scintillation crystal 125 and 
photomultiplier tube 127 are axially spaced at a distance 131 from the source 120. Neufron and 
gamma ray shielding 116 is disposed between the long and short spaced detectors. As in the 
borehole embodiment shown in Fig. 1, it is preferred that the first and second scintillation 
detectors be identical in response characteristics to facilitate specfral normalization. Both the short 
and long spaced detectors are displaced radially from the major axis of the tool 10'. 
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The tool 10' is positioned against or near the material 1 14 to be measured as shown in Fig. 
7. Between the long and short spaced detectors and the material 1 14 to be measured, the gamma 
ray and neutron shielding material 116 is fabricated to form angled collimation conduits defined 
by the walls 119. The collimator conduits are filled with thermal neutron shielding material 122. 
This enhances sensitivity to Compton scatter radiation by coUimating the scintillators to preferred 
angles of scatter. Furthermore, the radial displacement of the two source-scintillator pairs away 
fi-om the material 1 14 being measured allows the thickness of material between the pairs and the 
material to be maximized for a given radial dimension of the tool 10'. 

Data collection methods, processing methods and auxiliary equipment for the system embodied as 
shown in Fig. 7 are the same as those for the system embodied as shown in Fig. 1 and in Fig. 2. 

To summarize, the bulk density measuring system utilizing neutron-induced high energy 
gamma radiation which is dispersed within the measured material provides an increased depth of 
investigation when compared with prior art systems using lower energy isotopic gamma-ray 
sources localized to the tool. This increased depth of investigation is obtained without 
necessitating the use of an accelerator source, which is costly to fabricate, costly to maintain, and 
operationally unreliable when compared to "chemical" or isotopic neutron sources. In addition, 
the system generates inelastic scatter, thermal neutron capture, and elemental activation gamma 
radiation which can be used to determine a variety of other material properties such as elemental 
concentrations, lithology, porosity and the like. 

While the foregoing disclosure is directed toward the preferred embodiments of the 
invention, the scope of the invention is defined by the claims, which follow. 
What is claimed is: 
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